*------------------------------------------------------------------------------*
* Family history - Dinas, Fouka, Schläpfer - JoP
*------------------------------------------------------------------------------*

*------------------------------------------------------------------------------*
* Figure D7
*------------------------------------------------------------------------------*


*------------------------------------------------------------------------------*
* Greece
*------------------------------------------------------------------------------*
use greece, clear

gen inter=asiaminor*T
local controls female i.age i.edu i.income i.occupation i.nomos 

estimates clear
eststo proud: reg proud_binst T `controls' if asiaminor==0, cl(municipality)
eststo id: reg id_nationalityst T `controls' if asiaminor==0, cl(municipality)

coefplot (proud, offset(-0.01) aseq("Proud to be Greek") ciopts(lcolor(blue blue)) mcolor(blue)), ///
levels(95 90) xline(0, lcolor(gs10)) keep(T) ///
ysca(titlegap(2)) msize(small) title("Greece", size(small)) ///
xtitle("", size(small)) xtitle("Treatment effect", size(small)) ///
xlabel(, nogrid labsize(small)) ylabel(, nogrid labsize(small)) xsca(titlegap(2)) ///
legend(region(lcolor(white))) swapnames legend(off) grid(none) 

*(id, label(Identity: nationality) ciopts(lcolor(blue blue)) mcolor(blue)), ///

graph save Greece, replace

*------------------------------------------------------------------------------*
* Germany
*------------------------------------------------------------------------------*

use germany, clear

local controls i.sex i.state i.age i.F4 i.F13
estimates clear
eststo proud: reg proud_german_st t1 `controls' if expellee==0, cl(GKZ)

coefplot (proud, offset(-0.01) aseq("Proud to be German") ciopts(lcolor(blue blue)) mcolor(blue)), ///
levels(95 90) xline(0, lcolor(gs10)) keep(t1) ///
ysca(titlegap(2)) msize(small) title("Germany", size(small)) ///
ytitle("", size(small)) xtitle("Treatment effect", size(small)) ///
xlabel(, nogrid labsize(small)) ylabel(, nogrid labsize(small)) xsca(titlegap(2)) ///
legend(region(lcolor(white))) swapnames legend(off) grid(none) 

graph save Germany, replace


*------------------------------------------------------------------------------*
* Combine
*------------------------------------------------------------------------------*

graph combine Greece.gph Germany.gph, xcommon iscale(2) ysize(2)
